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ABSTRACT 



Through the analysis of a set of numerical simulations of major mergers between initially non-rotating, pressure supported progenitor 
galaxies with a range of central mass concentrations, we have shown that: (1) it is possible to generate elliptical-like galaxies, with 
v/cr > 1 outside one effective radius, as a result of the conversion of orbital- into internal-angular momentum; (2) the outer regions 
acquire part of the angular momentum first; (3) both the baryonic and the dark matter components of the remnant galaxy acquire part 
of the angular momentum, the relative fractions depend on the initial concentration of the merging galaxies. For this conversion to 
occur the initial baryonic component must be sufficiently dense and/or the encounter should take place on a orbit with high angular 
momentum. Systems with these hybrid properties have been recently observed through a combination of stellar absorption lines and 
planetary nebulae for kinematic studies of early-type galaxies. Our results are in qualitative agreement with such observations and 
demonstrate that even mergers composed of non-rotating, pressure-supported progenitor galaxies can produce early-type galaxies 
with significant rotation at large radii. 
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1. Introduction 

It has been known for several decades that early type galax- 
ies have complex kinematics and varying amounts of rotation 
and that these c haracteristics depend on lum inosity and isopho- 
tal shape (e.g.. lDavieslll983l: lBendeij|l98 8'). Rece ntly, the re- 
sults of the SAURON survey jBacon et al. 12 001 ; de Zeeuw et all 
I2002h have lead to a division of early-type galaxies into two dis- 
tinct classes, slow and fast rotators, depending on the amount 
of the angular momentum per unit ma ss of their stellar com - 
ponent inside one effective radius (Re-, 'Emse llem et alj|2007[) . 
These two classes also have different properties, slow rotators 
are g enerally more massive , less flattened systems than fast rota- 
tors jCappellari et al.ll2007l) . Studies of the more distant regions 
of early-type galaxies, traced by planetary nebulae, have shown 
that their complexity of the kinematics of the inner regions also 
extends to the outer halos (Co ccato e t al. 2009). Interestingly, 
these outer halo es can be more rotat ionally dominated than their 
central regions dCoccato et al.ll2009h . 

Extensive studies using numerical simulations suggest that 
different processes lead to the formation of slow and fast 
rotating early-type galaxies. Most of these studies have fo- 
cused on the remnants of spiral-spiral mergers, which is 
likely to be an important process for the evolution of forma- 
tion of ellipticals and perhaps responsible for driving galax- 
ies from the blue cloud to the red seque nce dCattaneo et alj 
120061: iFaber et al.ll2007l: iRomeo et al.ll2008h . Even if the agree- 
ment between simulations and observations is often qualita- 
tive only - which is most likely due to incomplete knowledge 
of the characteristics of the progenitors of these merge r rem- 
nants - sorne con c lusions can still be d rawn ( Ouinn et a lJI 19931: 
iNaab et all [19991: IVelazquez & White! [T999i: ,Bendo & BarnesI 



2000';' Cretton et al .'200l';'Naab & Burkert'2003VBournaud et al 
2004, 2005; Naab et al.. .2006a; Gonzalez-Gai-cfa & BalcelTs 
2005; Gonzalez-Garcia.etalJ 120061; iJesseit et al.1 l2007l |2008|) 



Mergers of two spirals can produce rotationally supported early- 
type galaxies (at radii inside 1 Rg), depending on the mass ra- 
tio of the progenitors. Major mergers (1:1 mass ratio) are suf- 
ficiently violent to completely destroy the initial ordered ro- 
tational motions in the progenitor disk£], typically producing 
a pressure-supported system. Apparently, mergers with higher 
mass ratios can generally produce remnants with significant ro- 
tation. Virtually all the focus so far has been on simulated spiral- 
spiral interactions, attempting to answer to the question: Is it 
possible to (at least partially) preserve the initial internal angular 
momentum of the progenitor disks? 

Encounters between pressure supported (spheroid dom- 
inated ) galaxie s have been systematically studie d by Whitj 



19781 7l979h; 



'Gonzalez -Garcfa & van Albadal (l2005albt) ; 

iGonzalez-Garcfa et a l. {200^. These studies show that such 
mergers r esult in a wide variety of remnant p roperties. In 
particular, iGonzalez-Garcfa & van Albadal (12005 al) pointed out 
that mergers of spherically symmetric galaxies, without dark 
matter, can result in remnants partially supported by rotation 
and that part of the or bital angular momentum is absorbed b y 
the halo, when present ( iGonzalez-Garcfa & van Albadall2005bh . 
However, none of these studies have investigated in detail 
how the orbital angular momentum is redistributed during the 



' Unless peculiar confi gurations are considered l iPfennige3 1 19971 : 
iPuerari & Pfennigei 1200 ih or the amount of gas in the progenitor disks 



is sufficiently high i Springel & Hernauislll2005l ; [Robertson et al.ll2006l : 
iHopkins et aLll2009h . 
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interaction, nor have they investigated the dynamical properties 
of the remnants out to large distances (several Re). 

In this Letter, we show, for the first time, that mergers of ini- 
tially spherical, pressure-supported stellar systems without ro- 
tation can produce hybrid remnants having elliptical-like mor- 
phology, but are dominated by rotation (v/cr > 1). As we shall 
demonstrate, the final characteristics of the remnant depends on 
both the initial orbital angular momentum and the central density 
of the merging galaxies. 

2. Models and initial conditions 

We study the coalescence of two equal-mass elliptical galaxies, 
consisting of a stellar and a dark matter (DM) component, dis- 
tributed in a Plummer density profile, without any dissipational 
component and initial rotation. The density profile of the dark 
matter halo is the same in all the simulations, we only changed 
the central density of the baryonic component in order to study 
how changing the density affects the angular momentum (AM) 
redistribution during the encounter. We thus considered gEOl, 
gEO and gEOm models, with increasing central stellar density 
(Table [1] and Fig. [1]), whose half-mass radii rgo and effective 
densities cr^/rso' are in agree ment with ob servational estimates 
for galaxies with M,- - -22 ( iDesroches a l. 2007), assuming a 
M/L=3 in the r-band . Since the total mass of the stellar com- 
ponent is fixed, increasing the central stellar density results in 
a more concentrated and compact stellar profile. In these mod- 
els, the DM contributes ~ 10-30% of the total mass within the 
half-mass radius of the baryonic component, in agreement with 
observational estimates dBarnabe et al.ll2069h . The initial galax- 
ies have been then placed at a relative distance of 100 kpc, with 
a variety of relative velocities, in order to simu late different or- 
bits (o rbital energies and angular momenta; see'Di Mat teo et al] 
(|2009|) ). In particular, in the following we will show results from 
three of these orbits (id=01, 05, 15), whose main parameters 
are given in Table [2l When refering to specific encounters, the 
nomenclature adopted is the following: morphological type of 
the two galaxies in the interaction (gEO, gEOl or gEOm), plus the 
encounter identification string (see first column in Table 2). Each 
pair has been modeled with = 120000 particles, distributed 
among stars (A^^^^ - 80000) and dark matter (A^bm=40000). 

All the s imulations have been run u sing the Tree-SPH code 
described in lSemelin & CombesI ( |2002|) . A Plummer potential is 
used to soften the gravity at small scales, with constant softening 
lengths of e = 280 pc for all particles. The equations of motion 
are integrated using a leapfrog algorithm with a fixed time step of 
0.5 Myr. With these choices, the relative error in the conservation 
of the total energy is of the order of 2 x 10"^ per time step. Since 
the work presented here only investigates dry-mergers, only the 
part of the code evaluating the gravitational forces acting on the 
systems has been used. 

3. Results 

3.1. Angular momentum redistribution: dependence on the 
orbits 

The first question we would like to address is: how orbital AM 
is redistributed and converted into internal AM during an inter- 
action and how does this depend on the orbital parameters of the 
encounter? In Fig. |2] we show the evolution of the total, orbital 
and internal AM, for two different orbits (id=15 and id=05), cor- 
responding, respectively, to increasing initial orbital angular mo- 
menta. The AM is evaluated with respect to the barycenters of 



the two galaxies, considering all the particles. The main results 
are: 
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Fig. 1. Initial properties of the galaxy models: (a) projected den- 
sity map of the stellar component of the gEO galaxy; (b) volume- 
density profiles of the galaxy models. The model is as indicated 
in the legend of the figure. 



Table 1. Parameters for the standard (gEO), less concentrated 
(gEOl) and more concentrated (gEOm) model galaxies. Both the 
stellar and dark matter profiles are represented by Plummer mod- 
els, having characteristic masses, respectively, given by Mg and 
Mh for the baryons and halo respectively, and core radii given by 
rB and for the baryons and halo respectively. The initial half- 
mass radii rgo of the baryonic component and effective densities 
cr jrsQ^ are also given. 







rB,rH 


'■so 






[2.3 X 10" Mq] 


[kpc] 


[kpc] 


[(km/s/kpc)2] 


aEO 


70, 30 


4.0, 7.0 


5.2 


3.2 


gEOl 


70, 30 


6.0, 7.0 


7.5 


2.9 


gEOm 


70, 30 


2.0, 7.0 


2.8 


3.9 



Table 2. Orbital parameters 



orbit id 


[kpc] 


[100ytm.r'] 


[Wkms-^kpc] 


[Wkm^s-^] 


01 


100. 


2.0 


57.0 


0. 


05 


100. 


2.0 


80.0^ 


0. 


15 


70.0 


1.6 


41.3 


-1.57 



" The initial distance between the two galaxies. 
The absolute value of the initial relative velocity between the two 
galaxies. 

The absolute value of the orbital angular momentum of the unit 
mass, i.e., L =\ Viai x Vj„i |. 
The total energy of the relative motion, i.e., 
E = Vi,/I2 - G(mi + in2)/ri„i, with mi = m2 = 2.3 x 10"Mo. 
° Even if the initial relative velocities in orbit 05 and 01 are equal, 
their orbital angular momenta differ because the initial orbital tan- 
gential velocities are different. 



- The total AM is (obviously) conserved during the interac- 
tion. 

- The orbital AM is constant until the first pericenter passage 
between the two galaxies. When this happens, dynamical 
friction and tidal torques act efficiently on the systems, con- 
verting part of the orbital AM into internal AM. 
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- At every successive close passage between the two galaxies, 
part of the orbital AM is converted into internal rotation of 
the two systems. This process ends when all the orbital AM 
is converted into internal AM and the two galaxies finally 
merge. 

- Instead of defining the merging time as the point when the 
two galaxy centres are sufficiently close we suggest, based 
on dynamical arguments, that the merging time could be the 
point when orbital AM is totally converted into internal AM 
(as we verified, these two definitions give similar merging 
times). This definition, of course, applies only for orbits hav- 
ing an initial orbital AM different from zero. 

- Because of AM conservation, the final amount of internal 
AM of the remnant galaxy depends on the initial amount of 
orbital AM. 

In Fig. [3] we show the evolution with time of the specific 
internal AM for both the stellar and the dark matter components. 
For each of the two galaxies, and for each component (baryons 
and halo), we evaluate the specific internal AM in four different 
radial regions. From this analysis, we can deduce that: 

- The orbital AM is converted into internal AM "outside-in", 
the external regions acquire part of the angular momentum 
first. 

- Interestingly, between the first pericenter passage and the fi- 
nal phases of coalescence, the outer regions show some ro- 
tation, while the inner regions do not. 

- Both the baryonic and the dark matter component acquire 
part of the orbital AM, and both the components, in the final 
remnant, show some rotation. The two initially non-rotating 
and spherically symmetric components are thus transformed 
into a flattened rotating systen^ 

- The higher the initial orbital AM, the higher the internal spe- 
cific AM of the final remnant, at all radii. 



Time [Myrl 



Time [Myrl 



Fig. 2. Evolution of the total (black line), orbital (solid red line) 
and internal (dashed red line) AM for two different orbits. The 
two galaxies are identical, so the internal and orbital AM are 
shown only for one of the two. The angular momenta are all 
measured perpendicular to the orbital plane and are in units of 
2.3 X IQi'Mokpckms-'. 



3.2. Angular momentum redistribution: dependence on the 
galaxy central density 

To investigate how the evolution and redistribution of AM de- 
pend on the initial central density of the two progenitors, we 
also simulated mergers with a range of concentrations (gEOl- 
gEO - less concentrated and gEOm-gEOm - more concentrated). 



^ As we confirmed, the remnants show some velocity anisotropy and 
both rotation and anisotropy contribute to the flattening. A discussion 
on this will be presented in a more detailed paper. 




_gE0gE005 


1 






1 , 1 



1000 



2000 



3000 



Time [Myr] 



2000 
Time [Myr] 



4000 



Fig. 3. Evolution of the specific AM, I, for four different re- 
gions of one of the two elliptical galaxies (r < 2kpc, solid lines, 
2kpc < r < 5kpc dashed lines, 5kpc < r < IQkpc, dot-dashed 
lines, and IQkpc < r < IQkpc, dotted lines. The specific AM 
of the stellar component is shown in red, that of the dark matter 
component in gray. I is in units of 100 kpc km s In the two 
panels, the initial orbital parameters of the interaction has been 
varied, so to have an increasing orbital AM, going from left to 
right, while the morphology of the interacting galaxies has been 
kept fixed. The specific angular momenta are all measured per- 
pendicular to the orbital plane. 
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Fig. 4. Panels (a), (b) and (c) show results simular to Fig. [3] but 
now the panels show the evolution of the specific AM for inter- 
acting pairs having initially the same orbital parameters, but dif- 
ferent morphologies (the stellar central densities increase from 
panel a-c). Panel (d); Specific AM, /, of the remnant galaxies as 
a function of radii containing a fixed percentage of the baryonic 
(red lines) and halo (black lines) mass. / has been evaluated at 
least 1 Gyr after the coalescence of the two progenitors. Line- 
styles correspond to the different morphologies shown in panels 
(a), (b) and (c), i.e., gEOlgEOlOl (dot-dashed line), gEOgEOOl 
(dashed line), gEOmgEOmOl (solid line). 



Fig. |4] shows the evolution with time of the specific internal 
AM for four different regions in the galaxy (see also Fig. [3). 
Since the initial orbital AM are the same in these simulations, 
the amount of internal AM of the final remnant will also be the 
same. However, AM is redistributed in a quite different way, de- 
pending on the initial concentration of the progenitor galaxies. 
We find that mergers with higher initial galaxy central density 
produce remnants with larger relative amount of rotation at any 
radius within a radius of 20 kpc. This result is true for both the 
baryonic and dark matter components. 

This behavior is related to the fact that the three remnant 
galaxies have different mass distribution profiles, as we will dis- 
cuss in the next Section, and so ultimately the physical regions 
chosen in Fig. |4] contain different amount of baryonic and dark 
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matter mass. A more appropriate comparison is the distribution 
of the specific AM within a constant fraction of the mass. We 
show the distribution of the specific AM as a function of radii 
containing a fixed percentage of the baryonic (red lines) or dark 
matter (black lines) mass in Fig.|4] panel (d). Again, this is eval- 
uated at least 1 Gyr after the merger has completed. The internal 
rotation is distributed in a quite different way in the three rem- 
nants as a function of mass. The baryons in each of the models 
have a similar distribution of / inside the effective radius, while 
outside it, less concentrated progenitors lead to a higher amount 
of specific AM. The dark matter component, however, shows 
the opposite trend - the specific AM is higher for all values of 
the enclosed mass for progenitors with high concentration. This 
is because more concentrated baryonic components are less re- 
sponsive to tidal torques, while the outer, less concentrated halo 
is more strongly affected by the tidal interaction. To illustrate 
this, in the case of the gEOlgEOlOl interaction, at the end of the 
simulation about 67% of the orbital AM has been acquired by 
the baryons, and only 33% by the dark halo, while for the most 
concentrated case (gEOmgEOmOl) the two components have ac- 
quired similar percentage of the initial orbital AM (48% for the 
baryons and 52% for the halo). 

Thus baryonic components with low concentration are more sus- 
ceptible to tidal torques, and thus acquire a higher percentage 
of the initial orbital angular momentum. 

3.3. Mergers of dense pressure-supported systems: 
rotationally-dominated systems with an elliptical-like 
morphology 

The varying amount of specific AM found in remnants of pro- 
genitors having different initial concentrations is also reflected 
in the kinematics. A detailed analysis of the relationship be- 
tween the morphological and dynamical properties of mergers 
of pressure-supported systems will be the topic of a future paper. 
Here we simply note that when a sufficient amount of specific 
AM is imparted to the central regions of the remnant galaxy, 
such remnants may show interesting hybrid properties. Some 
may have a morphology typical of an elliptical galaxy, but a v/cr 
ratio higher than 1, over most of their extended stellar distribu- 
tion. This is the case, for example, in the gEOmgEOmOl merger 
remnant. Even if its baryonic component has a specific AM in- 
side the half-mass radius similar to that of the less concentrated 
galaxies (panel d of Fig.|4]i, the gEOmgEOmOl remnant is much 
denser than, for example, the gEOlgEOlOl remnant (Fig.|5]l. The 
different mass distribution of the two galaxies is reflected in dif- 
ferent dynamical properties in an interesting way: even if they 
acquire a lower percentage of the orbital AM, being less sus- 
ceptible to tidal torques, the baryonic components of the most 
concentrated systems show a higher central velocity dispersion, 
a steeper dispersion profile, a higher value of the line-of-sight ve- 
locities, and a \'/cr ratio which becomes greater than 1 at about 
r = rso (the baryonic half-mass radius; Fig.|6]l. This may appear 
paradoxical, but this is due to the fact that the internal AM is 
redistributed over a much more compact distribution compared 
to low-co ncentration remnants. Most of the galaxies studied by 
ICoccato et al. (2009) show a v/c ratio not greater than 0.6, while 
a fraction of them have kinematics that become increasingly sup- 
ported by rotation in the outer parts, qualitatively in agreement 
with the hybrid merger remnants discussed here. This supports 
the idea that the morphological properties of a galaxy are not 
univocally related to the dynamical ones. Mergers can indeed 
produce hybrid systems, h aving a spiral-like morphology but 
an elliptical-like kinematics (iJog & Chitrel2002l : lBoumaud et al.l 




r [kpc] 




Fig. 5. Morphological characteristics of two remnants of 
pressure-supported, initially non-rotating galaxies. Top pan- 
els: (Left) Projected density map of the remnant of the 
gEOmgEOmOl encounter (progenitors with high concentrations); 
(right) Surface density profile of the baryonic (red color) and 
dark matter component (black color) of the remnant (solid line) 
and its progenitor galaxy (dashed line). The green arrows in- 
dicate the half mass radius (rso) of the baryonic component. 
Bottom panels: Same as top panels, but progenitors with rela- 
tively low concentrations (gEOlgEOlOl). 




Fig. 6. Kinematical properties of two remnants of E-E merg- 
ers. Left panel: line-of-sight velocities (blue) and velocity dis- 
persions (black) of the gEOmgEOmOl (thick lines) and of the 
gEOlgEOlOl (thin lines) remnants. Right panel: corresponding 
v/o" ratio (gEOmgEOmOl, thick line; gEOlgEOlOl; thin line). 



I2004l2005|) or, as presented here, elliptical-like morphology but 
are rotationally-dominated. 

4. Conclusions 

Through the analysis of a set of numerical simulations of major 
mergers between initially non-rotating, pressure supported pro- 
genitor galaxies, we have shown that it is possible to generate 
elliptical-like galaxies, with v/cr > 1 outside one effective ra- 
dius, simply as a result of the conversion of the orbital AM into 
internal one, during the mergers of two initially non-rotating pro- 
genitors. For this to occur the initial baryonic component must 
be sufficiently dense and/or the encounter should take place on a 
high AM orbit. 

Systems with suc h hybrid prope r ties ha ve been recently ob- 
served. For example, ICoccato et al.l (|2009|) found that galaxies 
in their sample with the most steeply declining velocity pro- 
files have generally fast rotating halos, as found for our fast ro- 
tating remnants. Our results are in qualitative agreement with 
such observations and demonstrate that even mergers composed 
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of non-rotating, pressure-supported progenitor galaxies can pro- 
duce early-type galaxies with significant rotation at large radii. 
In particular, we stress the important point that this mechanism 
does not require initially any internal AM in the colliding galax- 
ies. 

Simulations by iNaab et all ( l2006bh show that repeated bi- 
nary early-type mergers lead to the formation of anisotropic, 
slowly rotating elliptical galaxies (inside R^). In this context, 
whether the hybrid systems presented in this Letter can be 
formed frequently and maintained also after successive mergers 
should be studied with a larger set of simulations, taking into ac- 
count repeated mergers (major and minor) with different relative 
inclinations, as well as progenitors with different DM profiles 
(lNavarroetal.1119971) . 

Nevertheless, it is possible to hypothesize a scenario where 
these elliptical, rotationally-supported systems are mainly found 
among remnants of dense early-type galaxies. These hybrid fea- 
tures will tend to vanish if successive major mergers take place, 
with a random orientation of the internal and orbital AM with 
respect to the spin of the "hybrid" galaxy. But this has yet to be 
investigated numerically. 
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